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Optically-trapped two-component Fermi gases near a broad Feshbach resonance exhibit universal 
thermodynamics, where the properties of the gas are independent of the details of the two-body 
scattering interactions. We present a global proof that such a universal gas obeys the virial theorem 
for any trapping potential U and any spin mixture, without assuming either the local density 
approximation or harmonic confinement. The total energy of the gas is given in scale invariant form 
by E = (edU/de), where e is an arbitrary energy scale in terms of which all length and energy scales 
that appear in the confining potential are written. This result enables model-independent energy 
measurement in traps that are anharmonic as well as anisotropic by observing only the cloud profile, 
and provides a consistency check for many-body calculations in the universal regime. 

PACS numbers: 03.75.Ss, 32.80.Pj 
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Strongly-interacting Fermi systems, such as neutron 
matter and resonantly interacting atomic Fermi gases, 
can exhibit universal behavior. These systems satisfy 
the unitary condition, where the zero energy scattering 
length a greatly exceeds the interparticle spacing, while 
the range of the scattering potential is neglig ible El US]. 
Unitary conditions are produced in an optically trapped 
Fermi gas by using a magnetic field to tune near 
a broad Feshbach resonance [a, IE 0j where strong in- 
teractions are observed According to the universal 
hypothesis, unitary systems must exhibit universal ther- 
modynamics 0, @, S H, S EH 11 1) since the interpar- 
ticle spacing sets the only natural length scale at zero 
temperature, so that the local properties can be written 
generally as functions of the density and temperature. 
As the properties of the system are independent of the 
details of the interactions, atomic gases can be used to 
test predictions for unitary systems in fields well out- 
side atomic physics. Recent model-independent studies 
of the entropy versus energy of a unitary Fermi gas 11211 
have been used to demonstrate universal behavior [13j |. 
Model-independent energy measurement is based on the 
virial theorem, which holds for a trapped unitary Fermi 
gas 0. 

We have demonstrated both theoretically and exper- 
imentally that the virial theorem holds for a trapped 
Fermi gas in the universal regime at a broad Feshbach 
resonance The proof is based on the local density 

approximation (LDA) with a local pressure P — 2£/3, 
which holds for a universal gas Here, £ is the local 
energy density (kinetic and interaction energies). Using 
the balance between the trapping and pressure forces, 
one easily obtains the total energy per particle E, which 
comprises the kinetic energy, interaction energy, and trap 
potential energy, 

E = (u+^-VUy (1) 

where U is the trapping potential. A scalar pressure re- 
quires that (xdU/dx) — (ydU/dy) — (zdU/dz). In an 



anisotropic harmonic trap, this leads to equal energies 
in all directions and E — imu>1{z 2 ) where m is the 
atom mass and lo z is the harmonic oscillator frequency for 
the z direction, which is precisely measurable. In prac- 
tice, the long axial direction of the cloud is most easily 



imaged 12j. Within the LDA, Eq. Q] is readily modi- 



fied to include corrections arising from anharmonicity in 
the trapping potential (l2l |. Amazingly, at resonance, a 
complex many-body Fermi gas, generally containing con- 
densed (superfluid) atom pairs, noncondensed pairs, and 
unpaired atoms, all strongly interacting in the nonper- 
turbative regime, obeys the same virial theorem as an 
ideal, noninteracting gas. 

Unfortunately, the local density approximation can 
break down in several cases of importance, such as quan- 
tum confined systems. In standing wave optical traps, for 
example, the harmonic oscillator frequency for vibrations 
along the standing wave direction z can be so large that 
the axial energy scale Tiuj z is large compared to the en- 
ergy in the transverse x—y directions. In spin-imbalanced 
mixtures, phase separation occurs and a surface tension 
exists [15|,|l6|. At low temperatures, the cloud separates 
into a superfluid core containing equal densities of both 
spins and a region of excess noninteracting spins outside 
the core. The shape of the central core is affected by 
the unpaired spins, and the pressure is not necessarily a 
function of the temperature and local density, since the 
spin fractions vary rapidly in space. The axial dimension 
is therefore not necessarily equivalent to the transverse 
dimensions of the trap. To avoid the local density approx- 
imation, a global proof of the virial theorem is required. 

A global proof of the virial theorem for a universal 
gas has been provided by Werner and Castin (l7j |. who 
note that the proof also holds for anisotropic traps, a re- 
sult due to Chevy [l7j| . Recently, an extremely simple 
global proof of the virial theorem for universal gases has 
been given by D. T. Son, based on the finite temperature 
Hellman-Feyman theorem [3]. This proof assumes har- 
monic confinement, and is obviously valid for anisotropic 
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harmonic traps and for spin-imbalanced mixtures. How- 
ever, it does not allow for anharmonic trapping poten- 
tials. Tan has derived both the pressure and the virial 
theorem for a strongly interacting Fermi gas with a finite 
scattering length in a generalized confining potential [l9| . 
The results show that the correction is oc 1/a, as con- 
firmed by Braaten and Platter [2(| , who obtain the same 
result as Tan by using the Hcllman-Fcynman theorem. 
One expects that the virial theorem also will hold for a 
universal gas with three body resonances, as discussed 
recently by Nishida, Son, and Tan 21 1. 

We now provide a simple global proof of the virial the- 
orem in the universal regime valid for any spin mixture 
and any trapping potential. The result enables a general 
method for precisely measuring the total energy of a uni- 
versal gas in anharmonic trapping potentials. Our proof 
follows closely the method employed by Son [lij], and 
yields a general, global, scale-invariant relation between 
the energy and the trapping potential, without assuming 
harmonic confinement. This permits the determination 
of small corrections for anharmonicity in the trapping 
potential, which is essential for precision measurements 
of the energy. 

First, consider a Hamiltonian that is a function of a 
parameter A, i.e., H(X). Then, the expectation value 
of H(X) with the n th eigenstate yields (n, X\H(X)\n, A) = 
E n (X), where E n {\) is the eigenvalue. Taking the deriva- 
tive of this equation with respect to A and noting that the 
states remain normalized as A varies yields the Hellman- 
Feynman theorem for one eigenstate, dE n (X)/dX = 
(n,X\dH(X)/dX\n,X). 

We assume that the evaporation rate from a trap of 
finite depth is sufficiently slow that the system is in quasi 
equilibrium, so that the probability for each many-body 
quantum state is determined by a Boltzmann factor for a 
fixed total number of atoms N. Then, following Son [l8| . 
we consider the free energy, F = —ksTlnZ, where Z is 
the partition function, Z — S„ exp[— j3 E n (X)], and f3 = 
l/(fcflT). Differentiating F with respect to A and using 
the Hellman-Feynman result for each eigenstate easily 
yields the finite temperature Hellman-Feynman theorem 
for any parameter A, 



dF(X) _ I dH(X) \ 
\ dX J 



dX 



(2) 



The Hamiltonian H generally comprises the total ki- 
netic energy operator, which introduces the parameter 
ti 2 /m with m the atom mass, the interaction potential, 
which introduces the s-wave scattering length, and the 
trapping potential energy U = ^^trapi^i) with the 
position of the i th atom. In the universal regime, the de- 
tails of the resonant two-body interaction potential can- 
not appear in the properties of the gas: Compared to 
the interparticle spacing, the zero-energy s-wave scatter- 
ing length a for opposite spin atoms is infinite, and the 



range R of the potential is zero. The corresponding en- 
ergy scales are ft 2 /(ma 2 ) and h 2 /(mR 2 ), which are zero 
and infinity, respectively. For identical spin atoms, the 
scattering length is zero. Hence, for both resonantly in- 
teracting and noninteracting atoms, only the trapping 
potential U introduces natural energy scales. 

We choose one arbitrary energy scale = e, for exam- 
ple, 1 joule, and write all other energy scales ti that 
appear in the trapping potential in terms of e, i.e., 
6i = cti€, as described further below. Now, let A = e 
in Eq. and note that (dH/de) = (dU / de) , since only 
the trap potential can contain the energy scale e. Then, 
dF(e)/de = (dU/de). 

In the universal regime, the free energy F can be writ- 
ten in terms of the chosen energy scale, and it must be 
of the form F = e /(fc^T/e), where / is a dimensionless 
function of the temperature T . Note that we suppress 
all other variables, such as the total atom number N, 
the spin fractions, N\/N, the on etc, which are held con- 
stant. Differentiating F with respect to T, we obtain 
the entropy S = -dF/dT = -k B f(k B T/e), where /' 
denotes the derivative of / with respect to its argument. 
Then, edF/de = F-k B T f = F+TS = NE, where E is 
the total energy per particle. Hence, we obtain a simple, 
general, and global relation between the total energy per 
particle and the trapping potential, 



E 



dU 



(3) 



where the brackets denote a statistical and quantum av- 
erage with the density operator. Eq. [3] is obviously inde- 
pendent of the energy scale e. 

As a simple example, consider a power law trap- 
ping potential of the form U n = Uo(r / at ra p) n , where 
r is the position of an atom of mass m. To obtain 
normalizable states, we assume a confining potential, 
where Uo > for n > and Uo < for n < 0. 
Let Uo — ct\t and h 2 /{ma 2 rap ) = o^e. Then, U„ = 

e 1+n ^ ai a^ /2 (mr 2 /h 2 ) n / 2 . Eq. [3] immediately yields, 

E=(n/2 + l)(U n ), (4) 

which is the standard result. 

Now consider an isotropic gaussian trap, for which the 
single particle trap potential energy is Ut ra p — Uq[1 — 
cxp(—r 2 /a 2 rap )}. Writing the energy scales associated 
with the trap depth and 1/e width as above, we have, 



Utrap — Ot\E 



1 — exp 



-a 2 e ■ 



For small displacements, the trap potential is har- 



monic: Utrap = 



2 /2 with muj 2 



2Uo/a 2 rap . Eq. [3] immediately yields the virial theorem 
mCu 2 r 2 }. Note that the oscillator energy 



trap 

result, E/N 
is htu = 
scales. 



v2AiA2e and therefore introduces no new energy 
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Applying Eq. [3] to the gaussian potential, and keep- 
ing up to lowest order in the anharmonic corrections, we 
obtain the energy per particle 



, _ 2 2 3 (mwV)' 
hj = ( muj r 



8 U n 



Here the averages are determined experimentally by mea- 
suring the cloud profiles for both spin states in spin- 
imbalanced mixtures. This reproduces our previous ex- 
pression for the anharmonic correction 12j, which was 
obtained assuming the local density approximation with 
a scalar pressure. For an anisotropic gaussian trapping 
potential, we simply replace u 2 r 2 by u 2 x x 2 + Lu 2 y 2 + u 2 z 2 
and use u) = (uJxUjyUJz) 1 ^ . The proof goes through un- 
changed. 

For a nearly harmonic trap with cylindrical symmetry 
in x and y, the energy is readily determined from the two 
dimensional x — z cloud profile. The energy per particle 



E = 2mJ 2 x (x 2 ) + muj 2 z (z 2 }. 



(5) 



Here we do not assume that the axial and transverse 
contributions are identical, as required in the local den- 
sity approximation, since the local density approximation 
may not hold. 

In general, the quartic terms for the trapping potential 
can be measured. This is accomplished by determining 
the parametric resonance frequencies and the trap aver- 



age values of (z 



and (z 2 x 2 ) from the cloud profiles 



for an ideal Fermi gas tuned near the zero crossing where 
the scattering length is zero, i.e., near 528 G for 6 Li. For 
a trap with near cylindrical symmetry, we expect that 
the total potential has symmetry under uj x x <-> u) v y. For 
a given atom number, the energy is varied by release and 
recapture as described previously [12|. The dependence 
of the frequency on the quartic averages of the spatial co- 
ordinates, which are determined from the cloud profiles, 
will provide the appropriate anharmonic corrections to 
Eq.[U 
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